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a : 

We describe the decay out of a superdeformed band using the methods of reaction theory. 
' Assuming that decay-out occurs due to equal coupling (on average) to a sea of equivalent 

t-H , chaotic normally deformed (ND) states, we calculate the average intraband decay intensity 

and show that it can be written as an "optical" background term plus a fluctuation term, in 
C^l ' total analogy with average nuclear cross sections. We also calculate the variance in closed 

form. We investigate how these objects are modified when the decay to the ND states occurs 
ON . via an ND doorway and the ND states' statistical properties are changed from chaotic to 

regular. We show that the average decay intensity depends on two dimensionless variables 
in the first case while in the second case, four variables enter the picture. 

(N 1 

t-H , 

■ 81. Introduction 

O; 

| The first superdeformed (SD) rotational band to be observed was that identified 

in the nucleus gg 2 Dygg by Twin et air' A sequence of nineteen 7-rays of nearly 
constant spacing was observed which permitted the attribution of the moment of 
inertia of a symmetric prolate rigid rotor with elliptic axes in the proportion 1:1:2. 
Since then 320 SD bands have been observed in various mass regions extending 
from vl~20 to yl~240. 2 ) The stability of nuclei at any deformation can be related 
to the existence of energy gaps between shells of independent particle states. The 
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shell gaps which appear for certain nuclear deformations give origin to non-spherical 
configurations of special stability. 3 )' 4 ) 

As explained by Lopez-Martens et al., 5 ^ the 7-spectrum of a compound nucleus 
whose decay path includes an SD band typically contains 7-rays which result from 
four distinct stages: 

I Statistical 7-rays which populate the SD band after formation of the nucleus 
under consideration by particle evaporation following a fusion evaporation or 
fragmentation reaction. 
II Collective E2 7-rays from decay along the SD band. 

III Statistical 7-rays from excited normally deformed (ND) configurations. 

IV Collective E2 7-rays from decay along the yrast ND band. 

In cases where this kind of decay scheme applies, stages I and III correspond to 
the cooling of a hot and chaotic system while stages II and IV correspond to the 
decay of a cold and regular system. Other decay schemes which include SD bands are 
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possible. An SD state may fission or emit particles instead of decaying to the ground 
state by 7-emission. Stage II shows the following intensity profile as a function of 
transition energy, rotational frequency or spin: 6 ) gradual population at high spin, a 
plateau at intermediate spin and sudden depopulation at low spin (between 10 and 
30 h depending on the mass region). 

SD bands by definition are found in the second minimum of the potential energy 
surface in deformation space (hyperdeformed bands in the third if they can be ob- 
served). For a certain interval in spin this minimum coexists with the first minimum 
where the ND states are found. Between these two minima there is a barrier of 
typically a few MeV in height. Rotational bands exist in both minima over a large 
interval in spin. For low spins an ND band is yrast while for higher spins an SD band 
is yrast. At a certain spin these two bands cross. Decreasing in spin from the crossing 
point the excitation energy of the SD band relative to the ND band increases due to 
the different moments of inertia of the two bands. At some moment states in the SD 
band begin to decay to states in the ND minimum. This transition between stages II 
and III proceeds by tunneling through the multidimensional barrier in deformation 
space. 7 )' 8 ) 

In most cases it has not been possible to connect SD bands to the rest of 
the known decay scheme and consequently only 7-transition energies, lifetimes and 
branching rations are known. The reason is the complexity of the ND decay scheme 
and our incomplete knowledge of it. In just 10% of the 320 bands identified have 
absolute excitation energies, spins and parities been designated. 9 )~ n ) 

In the present paper we are interested in stage II of the 7-intensity. In fact 
we shall not discuss the feeding 12 ) of SD bands and restrict ourselves to the decay 
out. In Section [21 we present analytic formulae for the energy average (including 
the energy average of the fluctuation contribution) and variance of the intraband 
decay intensity of a SD band in terms of variables which usefully describe the decay- 
out. 13 )~ 15 ) In agreement with Gu and Weidenmiiller 14 ) (GW) we find that average of 
the total intraband decay intensity can be written as a function of the dimensionless 
variables r^/Ts and / D where i - ^ is the spreading width for the mixing of an 
SD state with ND states of the same spin, D is the mean level spacing of the latter 
and Ig (.T/v) are the electromagnetic decay widths of the SD (ND) states. Our 
formula for the variance of the total intraband decay intensity, in addition to the 
two dimensionless variables just mentioned, depends on the dimensionless variable 

r N /(r s + r L ). 

The results of Section |2] depend on two statistical assumptions. Firstly it is 
assumed that the states in the ND minimum are chaotic and that their statistical 
properties can be described by the GOE. Secondly, it is assumed that on average the 
SD band couples with equal strength to all ND states. However it might be appro- 
priate in some mass regions where the decay out occurs at lower excitation energy to 
describe the ND states by Poisson statistics. Further, it may be incorrect to assume 
that all ND states are equivalent. Certain ND states which have a stronger overlap 
with the SD states may serve as doorways 16 ) to the remaining ND configurations. 
In Section 0a model for which the ND states obey Poisson statistics and for which 
an SD state couples to a single ND doorway is discussed. 17 )' 18 ) It is shown that the 
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addition variables become important. 

§2. Energy average and variance of the decay intensity 

The total intraband decay intensity has the form 14 )' 17 ) 

/oo 
dE\A 00 (E)\ 2 , (2-1) 
-oo 

where Aqq(E) is the intraband decay amplitude and is the electromagnetic decay 
width of SD state |0). The energy average of Eq. Q2TJI may be written as the 
incoherent sum 15 )' 19 ) / = J av + Ifl, where 

"OO 



4v = 4v = (2vrr 5 )" 1 / d£ ^oo(^) (2-2) 

J — oo 



and 



/OO 
d£|Afl (ii:)| 2 , (2-3) 
-oo 

where we have written the decay amplitude as Aoq = Aqo + Aq where ^4oo is a 
background contribution and Aq is the fluctuation on that background. Sargeant 
et aZ. 15 ) took the background to be 

^ = E-Eo + i(T s + rl)/2- (2 ' 4) 

Eq. (|2-4j) exhibits the structure of an isolated doorway resonance. The doorway 
|0) has an escape width Jg for decay to the SD state with next lower spin and a 
spreading width r* = 2-kv 2 / D where v 2 is the mean square coupling of [0} to the 
ND states whose level spacing is D. 

The auto-correlation function of the decay amplitude is given by 15 ) 

A% (E)A* (E>Y « 2 (2^/D)- 1 (^/r s ) 2 A^E) 2 ^ _ .^ A^Effa) 



When E' = E this reduces to 



\AL\ = 2{2nr N /D)- 1 2 (2-6) 

which is the average of the fluctuation contribution to the transition intensity. 

The integrals in Eqs. ()2-2j) and (|2-3[) may be carried out using the calculus of 
residues. One obtains 

4 OE = (i + rVr 5 )- 1 , (2-7) 

for the average background contribution and 

h = 2 {<nr N /D)- 1 ( ri// ^) = 2 (irrN/Dy 1 I m (1 - / av ) 2 , (2-8) 

(i + n/r s y 
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for the average fluctuation contribution to the average decay intensity. Eq. (|2-8|) for 
Ifl was compared to the fit formula, obtained by GW, 



7 GW 



1 - 0.9139 (r N /D) 



0.2172 



exp < 



0.4343 In (£\ -0.45^) 



-0.1303 



2 



(r N /D)-°- 1477 



(2-9) 

and qualitative agreement of the two formulae found. 15 ) The dependence of Iq (and 
that of J av ) on r^/Ts results from the resonant doorway energy dependence of the 
decay amplitude Aqq(E) [Eq. I|2-4jl ] . This energy dependence also manifests itself 
in the average of the fluctuation contribution to the transition intensity |AQ (i?)| 2 
[Eq. (|2-6|) ]. GW include precisely the same energy dependence in their calculation 
by use of an energy dependent transmission coefficient to describe decay to the SD 
band. This is the reason for our qualitative agreement with GW concerning /. 

A measure of the dispersion of the calculated / is given by the variance (AI) 2 = 



(/ — J) 2 . It was shown 15 -* that 



(AT) 2 = J fl 2 /i (0 + 2J av Ifl/ 2 (f ) , (2-10) 
where the variable £ is defined by 

£ - = ^(i + r'/r*) = (2-ii) 

and 

MQ = WRi + (Thf + ^TW a " d Mi> = WTTy (242> 

Since the variance depends only on (Pg + r^)/r^ in addition to r^/Ts and 
r^/D, upon fixing the latter two variables the variance may be considered a function 
of any one of r l /T N , r s /r N , r^/D or r s /D [see Eq. CUD ]- Our result for the 

variance of the decay intensity, (AI) 2 [Eq. (|2-10jl ] has a structure reminiscent of 
Ericson's expression for the variance of the cross section. 20 ) In the case compound 
nucleus scattering, extraction of the correlation width from a measurement of cross 
section autocorrelation function permits the determination of the density of states 
of the compound nucleus. 21 ) In the present case the variance supplies a second 
"equation" besides that for I. Both equations are functions of and D, since the 
electromagnetic widths are measured. Thus both 1^ and D can be unambiguously 
determined. The derivation in this section is strictly valid only for r^/D S> 1. 15 ^ 
However the formulae were found to be qualitatively correct even for < l. 15 ) 

§3. Regularity versus chaos in the ND minimum 



Aberg 18 ) has suggested that an order-chaos transition in the ND states enhances 
the tunneling probability between the SD and ND wells and consequently that the 
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decay-out of SD bands is an example of "chaos assisted tunneling" . Crucial to the 
argument 17 ) is a postulate that the decay-out occurs via an ND doorway state. An 
ND doorway state can be visualized as being the tail in the ND minimum of the wave- 
function of the zero-point vibration in the SD minimum which may be constructed 
by the Generator Coordinate Method (GCM) of Bonche et al. 22 ^ Microscopically 
the vibration is a coherent superposition of lp-lh states which is damped by the 
two body residual interaction. Depending on the strength of the residual interaction 
the 2p-2h states will be regular, chaotic or intermediate between these two limits. 
Decreasing in spin from the point where the SD and ND bands cross, the excitation 
energy of the SD band relative to the ND band increases. Near the crossing point 
the residual interaction is weak, therefore the ND states may be characterized by 
quantum numbers which are approximately conserved and their energy spectrum ex- 
hibits degeneracies. In the language of quantum chaos the ND states are regular and 
obey Poisson statistics. As the spin decreases the residual interaction grows with the 
increasing excitation energy, destroying degeneracies and increasing the complexity 
of the ND states until the regime of quantum chaos is reached as characterized by 
the GOE. This is a plausible picture of the evolution of the nuclear Hamiltonian with 
decreasing spin. Given that it is correct, how is the cascade down the superdeformed 
band modified by an order-chaos transition in the ND minimum? 

Aberg 18 ) constructed a random matrix model for the ND states which inter- 
polates between Poisson and GOE statistics by varying the chaoticity parameter A 
continuously in the range < A < 1. The parameter A which is the determined 
by the ratio of the variances of the diagonal and off-diagonal random Hamiltonian 
matrix elements simulates the effect of the residual interaction. The SD state is 
assumed to lie in middle of the ND spectrum and the doorway is chosen to be an 
ND state which lies halfway between the SD state and the edge of the ND spectrum. 
The tunneling probability is given by Tsd^nd = Yl l°/J 2 where a'^ is the amplitude 
of the nuclear wavefunction in ND state fM. The broadening of the distribution of 
the a'^ by the residual interaction results in an increase in Tsd->nd which is denoted 
"chaos assisted tunneling" . The enhancement in the tunneling probability due to the 
transition from regularity to chaos is deduced from the ratio of Tsd-*nd for A = 0.001 
to Tsd^nd for A = 0.1. The value of this ratio is estimated to be N where N is the 
number of ND states. 

Although this idea is also plausible it is not clear what the effect of the change 
in Tsd^nd calculated by Aberg 18 ) has on the average of the total relative decay 
intensity of the collective E2 electromagnetic transitions of SD bands where the 
sudden attenuation is experimentally observed. We calculated 17 ' the background 
contribution to this quantity for Aberg's model using the reaction theory methods 
which were outlined in Section Inland conclude that the order-chaos transition cannot 
explain the decay-out. 

Instead of Eq. ()2-4|) let us now take the background decay amplitude to be 17 ) 



E-E + i[r s + W 00 (E)}/2' 
Following Aberg 18 - 1 we assume that |0) only mixes with one ND doorway state \d) 
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whose energy is Ed- The state \d) is subsequently mixed by the residual interaction 
with the remaining ND states, \Q), Q = 1, ...N, having the same spin as |0) and \d). 
The self-energy is then given by 

W 00 (E) = \V 0d ? f /r (3-2) 

^ Q E-E Q + ir N /2 

where Vo<i is the interaction energy of |0) and \d) and Cd(Q) is the amplitude of \d) 
in \Q). The \Q) lie in the interval L = ND where D here denotes the mean spacing 
in energy of the \Q). 

Ignoring an energy shift of |0), the background decay intensity becomes 



_ r s r JT7 , i 

'av — 



(3-3) 



2tt [E - £ ] 2 + [r s + 2^\VU 2 S d (E)f/^ 

where the doorway strength function is given by 

In Sargeant et al. 17 ^ the effect of the chaoticity of the ND states on I av was in- 
vestigated by varying the strength of the residual interaction of the \Q) and their 
interaction with \d), both being assumed to be proportional to the chaoticity parame- 
ter A. The limiting value A=0 results in the \Q) having Poisson statistics (regularity) 
while A=l results in their having GOE statistics (chaos). The value of A determines 
the shape of Sd(E) which is precisely the strength function that was previously 
investigated as a function of A. 23 ) 

Instead of studying the interpolation between the limits A=0 and A=l by nu- 
merically diagonalising random matrices and performing ensemble averages as we 
did before, 17 )' 23 - 1 here we restrict ourselves to the limiting case A=0. As A — > 0, 
c d(Q) $dQ so that Sd(E) reduces to the single Breit-Wigner term, 

sFd ° 1{E) = hjE^WTiyi- (3 ' 5) 

When A — > and i~jv — * 0, then Sd{E) — > 5(E — Ed). For non-zero A, Sd(E) broadens 
with increasing A until when A = 1 it takes a form, independent of /jy, which is well 
approximated by 23 ) 

qGOE, F \_f 1: \E — E d \ < L/2 , s 

Sd {E) -\ 0, \E-E d \>L/2 ■ (3 ' 6) 

Inserting Eq. (JSZEJ) into Eq. we find that Z^OE = (1 + r^/Ts)' 1 as long as 

in agreement with Eq. 1)2-7(1 of Section |2J Gu et al. discuss the 
strength function further for L finite 24 ) and L infinite. 25 ) Eq. (|3-3() for / av with 
Sd = Sj m depends on four parameters: Jgf, |Vbd| 2 , /jv and the distance in energy 
separating \d) from \0) , A = Ed — Eq. It is useful to introduce a spreading width 
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iog 10 ff7D) iog 10 (ryo) iog 10 (ryD) 

Fig. 1. The left graph shows the average background decay intensity versus log 10 (-Ts/D) [The 
solid line was calculated using Eq. 112 -711 and the dotted line using Eq. The middle graph 

shows the fluctuation contribution to the average decay intensity versus log 10 (Fs / D) [The solid 
line was calculated using Eq. 12- 811 with Eq. 12 -7> and the dotted line using Eq. 12-811 with 
Eq. lEZl- The dashed line was calculated using Eq. 12-91 1. The right graph shows the standard 

deviation of the decay intensity (=y (AI) 2 ) versus log 10 (.Ts/D) [The solid line was calculated 
using Eq. (12101) with Eqs. (12-811 and 12-7H and the dotted line using Eq. 12-1011 with Eqs. 12- 811 and 
ijPfjl . The dashed line was calculated using Eq. ll2~Tr])l with Eq. ((221] • The following numerical 
values were used to plot the graphs: r^/Ts = 1, Fn/D — 0.1 and A/D=l. 



now defined to be = 2ir\Vod\ 2 /D. Inserting Eq. (|3-5|) into Eq. Q3-3JI yields our 
result for J av in the Poisson limit: upon making the change of integration variable 
x = (E — Eq)/D this takes the form 

H- = ^Td X - ^ — (3-7) 



2tt 



It is seen from Eq. (H£ZJ) that, in addition to E l /E s , /f v oi is a function of a further 
three dimensionless variables: E^/D, E$/D and A/D. 

The introduction of the ND doorway is analogous to including another class 
of complexity in multistep compound preequilibrium reactions. 26 )' 27 -* Calculation 
of the fluctuation contribution to the average decay intensity and the variance in 
the presence of an ND doorway is beyond the scope of this paper. However, since 
Eq. (|2-8|) for the fluctuation contribution and Eq. (|2-1U|) for the variance are express- 
ible in terms of the background contribution, J av , one can obtain add hoc estimates 
by simply substituting Eq. (j3-7j) into Eqs. (|2-8|) and (|2-10|) (see Figure P). A rigorous 
calculation of the fluctuation contribution and variance when the ND states obey 
Poisson statistics may be possible using the averaging method of Gorin. 28 ) 
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The most important conclusion to be drawn from Figure^is that a measurement 
of the variance (standard deviation) of the decay intensity determines the density 
of ND states p = 1/D. For comparison we have also calculated the standard devi- 
ation by inserting the GW fit formula [Eq. (|2-9|l ] into our formula for the variance 
[Eq. (jlTDJ)]. 

Eq. (|3-7|) is formally identical to the branching ratio, Fs, calculated by Stafford 
and Barrett 29 ^ (SB). They model the decay-out by allowing a single SD state to mix 
with a single ND state. SB and Dzyublik and Utyuzh 30 ) who develop the method 
of SB further address the lack of knowledge of A by averaging over A assuming a 
uniform distribution. 31 ) Cardamone et a/. 32 ) on the other hand assume a Wigner 
distribution, Pqoe(^) = ttA/2 exp(— 7rZ\ 2 /4). The preceding discussion suggests 
that a Poisson distribution Pp Q i(Z\) = exp(— A) is more appropriate. 

An Addendum to Sargeant et al. 17 ^ is in preparation, where we shall discuss the 
question of whether the order-chaos transition can explain the decay-out when A is 
non-zero. 33 ^ 
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